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Abstract —In wireless fading channels, multi-user scheduling 
has the potential to boost the spectral efficiency hy exploiting 
diversity gains. In this regard, proportional fair (PF) scheduling 
provides a solution for increasing the users’ quality of experience 
hy finding a balance between system throughput maximization 
and user fairness. For this purpose, precise instantaneous channel 
state information (CSI) needs to be available at the transmitter 
side to perform rate adaptation and scheduling. However, in prac¬ 
tical setups, CSI is impaired by, e.g., channel estimation errors, 
quantization and feedback delays. Especially in centralized cloud 
based communication systems, where main parts of the lower 
layer processing is shifted to a central entity, high backhaul 
latency can cause substantial CSI imperfections, resulting in 
significant performance degradations. In this work robust rate 
adaptation as well as robust PF scheduling are presented, which 
account for CSI impairments. The proposed rate adaptation 
solution guarantees a fixed target outage probability, which is 
of interest for delay critical and data intensive applications, such 
as, video conference systems. In addition to CSI imperfections 
the proposed scheduler is able to account for delayed decoding 
acknowledgements from the receiver. 

I. Introduction 

T he steadily increasing throughput demand of mobile 
users together with the limited radio spectrum in cel¬ 
lular communication systems requires an efficient resource 
allocation. In orthogonal frequency division multiple access 
(OFDMA) where multiple users are served on different non¬ 
overlapping time-frequency resources, modern scheduling al¬ 
gorithms exploit multi-user diversity by accounting for the 
individual time variations of the users’ mobile radio channel. 

In order to perform scheduling, instantaneous channel infor¬ 
mation needs to be available at the transmitter side. However, 
channel state information (CSI) is typically impaired due to 
a variety of sources. In frequency division duplex (FDD) 
systems, CSI imperfections result from noisy pilot reception 
at the receiver and erroneous feedback transmission due to 
quantization errors as well as delays between channel obser¬ 
vation and transmission Cl, 0. The latter aspect is especially 
relevant at slow fading channels 0, where the provided CSI is 
still fairly correlated with the actual channel. Due to the lack of 
exact channel knowledge available at the transmitter, the rate 
achievable is not known precisely. This causes a non-optimal 
scheduling decision as well as imperfect rate assignment. If the 
channel amplitude is overestimated, decoding errors may occur 
at the receiver, which lead to outages of the respective data 
blocks and ultimately to a degradation in spectral efficiency 
0 . 

A. Related Literature 

In the area of multi-user scheduling, different objectives 
have been investigated over the last two decades. Opportunistic 



Fig. 1. Setup with K user equipments (UEs) served by a base station (BS), 
which is connected to a central node (CN). According to a cloud-based 
architecture, the main PHY and MAC layer processing is performed at a 
centralized data center. 


scheduling maximizes the system throughput by selecting 
the user with the best channel 0, 0. However, with this 
approach the system may be unfair towards users with poor 
channel conditions, which might never be scheduled. User 
fairness can be achieved by employing throughput balancing 
or round robin scheduling 0 which comes at the cost of 
a significant degradation in the overall throughput. A com¬ 
promise between throughput maximization and fairness is of 
interest O and has been reflected by the proportional fair 
(PF) scheduler, which maximizes the sum of the logarithmic 
user-throughputs 0, 0. PF multi-user scheduling has been 
investigated for several scenarios, such as, single-cell m and 
multi-cell networks im, ina. An extension to multi antenna 
systems has been discussed in m, m. A PF algorithm 
which accounts for intra-cell as well as inter-cell interference 
has been proposed by iia. 

The publications mentioned above inherently assume that 
CSI is precisely available at the scheduler. Scheduling with 
imperfect CSI has been studied in several works considering 
throughput maximization CD-lIISl. A PF scheduling scheme, 
which provides robustness against imperfect channel knowl¬ 
edge has been presented in our preliminary work ifT^ . 

For delay constraint applications, the system operates at 
a fixed outage probability and hybrid automatic repeat re¬ 
quest (HARQ) ||20l, im is subject to a maximum number 
of re-transmissions in order to guarantee a certain quality 
of service (QoS). For applications with strict conditions on 
latency requirements (e.g., video conference systems, tactile 
internet), HARQ may not be an appropriate mechanism. More¬ 
over, round trip delays may exceed respond time limitations, 
required by HARQ protocols (e.g. LTE-A). In LTE-A the 
transmission rate can be adapted by means of outer loop link 
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adaptation (OLLA), which is based on acknowledgments of 
previous transmissions. However, in scenarios where the round 
trip delay is in the range of the channels coherence time, 
also the acknowledgment is insufficiently correlated with the 
actual channel. In this work, we restrict our-selves to inner 
loop link adaptation (ILLA), where outage do not impact the 
rate adaptation of following transmissions. 

For generalized fading channels, outage probability has been 
studied in ll22l for receive combining techniques. Rate assign¬ 
ment for throughput maximization w.r.t. slow fading channels 
has been presented in ll23l . while fast fading channels have 
been considered in ll24l employing a fixed outage probability. 
The effect of imperfect CSI on the scheduling performance 
has been investigated in 1251. In 1251 . the authors presented 
an analytical framework for jointly analyzing the performance 
of certain state of the art methods regarding feedback transmis¬ 
sion, rate adaptation and scheduling. In contrast to our work 
the algorithms itself are not optimized. Instead the framework 
allows to select the optimal combination of particular methods 
for a given scenario. 

Recently, cloud based architectures in cellular communica¬ 
tion systems attract significant attention ll26l . The basic idea 
is the outsourcing of computationally complex algorithms of 
the physical (PHY) and the medium access control (MAC) 
layer to a data center which provides processing power as 
well as high scalability. In such an architecture, the scheduler 
is located at a central node (CN), which is accessible to a 
large number of BSs via the backhaul network. In practice, 
the backhaul connections cause technology-dependent delays. 
Such extra latency causes outdated CSI as well as delayed 
transmission acknowledgements. 

B. Contribution of this Work 

In this paper, we study the effect of imperfect CSI in 
the context of a cloud-based cellular communication sys¬ 
tem, where rate adaptation and scheduling is performed at a 
centrally located data center. In order to provide large user 
throughput together with a guaranteed QoS, a robust rate 
adaptation method is presented, which guarantees a fixed out¬ 
age probability. The improved rate adaptation is incorporated 
into the PF scheduling metric in order to achieve robustness 
against impaired channel knowledge. The algorithm addition¬ 
ally accounts for delayed acknowledgements of successful 
transmissions. 

This paper is an extension of our previous work lfT9l . It 
provides deeper insight on the derivations of proposed meth¬ 
ods regarding rate adaptation and PF scheduling. Moreover, 
algorithmic solutions w.r.t. a practical implementation as well 
as an extensive evaluation part is given. While the impact of 
backhaul delays has mainly been addressed in the context of 
cooperative cellular systems ||27l. ED, this work focuses on 
non-cooperative cloud-based communication systems. 

The remainder of this paper is structured as follows. The 
system model is presented in Section HIl before outage proba¬ 
bility results with imperfect CSI are given in Section |III] Sec¬ 
tion |IV] examines the proposed PF scheduling scheme, while 
Section |V] shows simulation results followed by conclusions 
in Section IVTl 
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Fig. 2. The message sequence chart illusttates the process between channel 
observation and transmission, including rate adaptation and scheduling at the 
central node (CN). 

C. Notation 

Conjugate, transposition and conjugate transposition are 
denoted by (•)*, and , respectively. Expectation is 
E{-}. The probability of an event A is P{A}, while the 
probability of event A conditioned on a given event B is 
P{A|i3}. Furthermore, C denotes the set of complex numbers 
and refers to a complex normal distribution with 

mean ^ and variance cr^. 

H. System Model 

This work considers a single-cell downlink system compris¬ 
ing K user equipments (UEs), which are served by the base 
station (BS) on different non-overlapping resource elements, as 
known from OEDMA. The BS is connected to a central node 
(CN), e.g., a data center with sufficient computational power 
and scalability. Accordingly, the main PHY and MAC layer 
algorithms are performed in a centralized fashion at the CN 
(see Pig. [III. In this work, the algorithms particularly comprise 
rate adaptation and multi-user scheduling. 

Centralized processing requires that instantaneous channel 
information is available at the CN. Assuming the employment 
of an PDD system, the downlink channel cannot be obtained 
by observing the uplink channel. Hence, CSI needs to be 
provided to the BS, from where it is then forwarded to the 
CN. The whole information flow including signaling as well 
as data transmission is illustrated as message sequence chart 
in Pig. [2] 

A. Data Transmission 

The instantaneous radio channel at time slot n between the 
BS and UP k is denoted with hk[n] ^ A/c(0, Afe), where 
Afc is the mean channel gain of UE k, reflecting path loss 
and shadow fading effects. A time slot reflects the duration 
of a transmission block, which refers to a collection of 
continuous symbol transmissions in time and frequency where 
the channel is assumed to be constant. The radio resources 
of each transmission block are scheduled only to a single 
user. Consequently, no inter-user interference needs to be 
considered. 











3 


Assuming Gaussian distributed channel input Xk ~ 
N'c{0,p) with transmit power p, the received signal is given 
as 

yk=hk[n]xk+rik, ( 1 ) 

where pk ~ A/c(0,1) is the zero mean complex Gaussian 
receiver noise at UE k with power equal to one. 

With precise knowledge of the instantaneous channel gain 
|/ifc[n]p at the transmitter as well as precise knowledge of the 
complex channel hk[n] (including amplitude and phase) at the 
receiver side, the capacity (maximum achievable transmission 
rate) for UE k at time slot n is given by 

Ck[n\ = log 2 (1 + p|/ifc[n]|^) . (2) 

Considering a long-term power constraint, the capacity is 
achieved by adapting the transmit power p at each time 
slot to the respective channel state 0. However, this work 
assumes an instantaneous power constraint with a fixed power 
allocation p at each time slot. The signal-to-noise-ratio (SNR) 
experienced at UE k is then given by jk = P^k- 

B. Outage Probability and Throughput 

The channel capacity in (|2]i is only achievable with precise 
rate allocation, i.e., the maximum rate supported by the 
instantaneous channel needs to be assigned for transmission. 
However, with imperfect channel knowledge, the achievable 
rate is not known at the transmitter. Consequently, the allocated 
transmission rate might not match the channel capacity. 

If the rate allocated is below capacity, the channel’s potential 
is not fully exploited. In this case only a reduced spectral 
efficiency is achieved. On the other hand, if the rate assigned 
exceeds capacity, outage of the transmitted data block occurs, 
i.e, the codeword received at the UE cannot be decoded 
correctly 0. 

Let function 5'fc[n] indicate the success of transmission to 
user k in time slot n, i.e., S'fc[n] is equal to one if the codeword 
has been successfully decoded and zero otherwise. With the 
rate assigned for transmission Rk[n], the actually experienced 
rate at UE k in time slot n is Rk[n] = Sk[n]Rk[n]. The 
probability of outage is given by 

Pout.fc W = P {Sk[n\ = 0} 

= IP {log 2 (1 + < Rk[n\} ■ 

Erom the transmitter point of view, hk[n] is a random variable. 
An extension of 0 to the case with side information is stated 
in Section Hn] 

The throughput of UE k obtained until time slot N is given 
by 

1 ^ 1 ^ 

Tk [^] = ;^ XI = ]v ^ 

n—1 n—1 

The throughput expression in (|4]i refers to the average rate of 
user k, experienced up to time slot N. 

C. Imperfect CSI 

In this section, we present the model for imperfect CSI used 
in this work. In EDD systems, the downlink channel cannot 


be obtained from corresponding uplink measurements but is 
observed at the UE side and provided to the BS via an uplink 
control channel with limited capacity. Consequently, the CSI 
available at the BS is impaired by channel estimation errors, 
feedback quantization and potential delays between channel 
estimation and actual data transmission. In this work, we use 
the model derived in ll2^ which considers all three effects. 
CSI is obtained by minimum mean square error (MMSE) 
estimation/prediction, i. e., the unknown channel 

hk [n] Me {hk [n] , (5) 

can be represented by a complex Gaussian random variable, 
whose mean value refers to the currently available channel 
estimate hk[n], while its variance Ck reflects the respective 
channel uncertainty 1^ . Note that Ck/M refers to the nor¬ 
malized CSI impairment and takes values in the interval [0; 1]. 
The unknown channel can equivalently be written as the sum 
of the channel estimate and a zero mean Gaussian random 
error ek[n] ~ A/c(0,efe), i.e. 

hk[n] = hk[n] + ek[n]. (6) 

Note that in this model, the mean gain of the channel estimate 
E{|/ifc[n]p} = Afe(l —Cfc/Afe) is reduced according to the error 
variance to Ck- According to ll29l . can be calculated based 
on the pilot SNR, the quantization resolution, the delay and 
the processing window length W of the prediction filter. In 
this regard the time variation of the channel is reflected by 
the correlation between two channel coefficients, delayed by 
A time slots: 

c[A] = En{hk[n]hl[n + A]}/Afc = Jo (^9^^ ^ 

where and Jq are the normalized coherence time (nor¬ 
malized to the time slot duration) and the zero-th order 
Bessel function of the first kind, respectively. The constant 
q = J|)"^(0.5) ensures that Tc is referring to the 50% coherence 
time. In this work, slow fading is assumed such that the 
channels of subsequent time slots show strong correlations. 
Using 0, the covariance vector 

c[A] = [c[A],...c[A + IU-l]]^ (8) 

as well as the covariance matrix C = [c[0],..., c[—A]] are 
defined, considering that W subsequent channel observations 
are used for prediction. 

Denoting the number of pilot signals within a transmission 
block as Np and the number of feedback bits per complex 
channel coefficient as Q, the error variance is given by 

cfc = (1 - 2-«)c^[A](C + pNp)-^c[A]. (9) 

A more detailed derivation and explanation of the feedback 
model is given in 1291 . 

III. Robust Rate Adaptation with Imperfect CSI 

In this section, the proposed rate adaptation scheme is 
derived. The algorithm is performed for each user individually. 
Since this work focuses on ILEA, the rate adaptation at a 
particular time-slot is not affected by previous transmissions. 
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Independence between multiple users and subsequent time- 
slots allows us to omit the UE index k as well as the time slot 
index n for improving readability in this section. Note, that the 
impact of feedback latency is still fully reflected by the error 
variance e. Furthermore, the amplitude of the actual channel 
and the CSI are denoted as g = \h\ and g = \h\, respectively. 

Utilizing the results of a, ED, the probability of outage 
conditioned on the available channel information g is given as 

Pout = p {log2 (l -b pg'^) < R \ g} 

= p{ 5 <^( 2 «- 1 ) /p I 3 } ( 10 ) 

= (v/( 2 ^-l)/p)t 

where Rg^^{x) denotes the cumulative distribution function 
(CDF) of g conditioned on side information g. Eq. (fTOl i 
allows to obtain the outage probability for the given channel 
estimate g when rate R is allocated. However, in case of delay 
constrained applications, the outage probability is given as a 
constraint and the transmitter is interested in the rate allocated 
to satisfy this constrained. Consequently, the inverse function 
Rig,Pout) needs to be found. 

The amplitude of a complex Gaussian non-zero mean ran¬ 
dom variable follows a Rician distribution 0 . Consequently, 
the probability density function (pdf) of g conditioned on g 
results in 

/»,*) = (-^) J„ (M) , (ID 

where Jo{x) = + 1)) refers to the 

modified Bessel function of the first kind and order zero, while 
r(a;) = /p is the gamma function. 

The CDF in (fTOl i is obtained by integrating (fTTl i over the 
interval [0; 6] with b = ( 2 -^ — l) /p: 

J^9\9ib)= t fg\gi9)d9- (12) 

Jo 

With the results of EH, the integral in (fT 2 T i is given by 
=l- 0 i 

= 1 _ exp (-£^) e::=o {(fjm {^), 

(13) 

where Qi is the Marcum Q-function and Jmix) = 
Sto z!r(i-4n-n) is the modified Bessel function of 

the hrst kind and order m. 

Hence, the resulting outage probability with given g is given 
by 

Pout = 1 -exp 

/ - \ ™ / A—\ 

V^oo ( g\/p \ r I 2gy/2^-l \ 

However, (O cannot be reformulated in order to determine 
the outage rate R in closed form. Hence, numerical methods 
need to be applied to determine the R which guarantees a 
target outage probability pout- 



Fig. 3. Rate assigned for transmission as a function of the amplitude of the 
available channel estimate. 



Fig. 4. Resulting outage probability as a function of the amplitude of the 
available channel estimate. 

Fig. |3] shows the rate allocated as a function of the estimated 
channel amplitude g. The blue curves refer to the proposed 
robust rate allocation (RA) scheme, which is given by the 
inverse function of (fT4l) . as discussed above. The graphs are 
plotted for a target outage probability of pout = 0.1 and 
Pout = 0.01. For comparison, non-robust RA is shown (black 
curve). The latter algorithm does not take into account that 
CSI is imperfect. Note, that with decreasing error variance e, 
the robust RA converges to the non-robust RA. The non-robust 
scheme can be parametrized by a back-off factor a, according 
to 

R = a\og2{l + pg^). (15) 

Decreasing a leads to a less aggressive rate allocation and 
a reduction of the outage probability. The results of the 
respective functions employing a back-off factor of a = 0.95 
are illustrated by the gray dashed line. 
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Fig. 5. Mean user throughput as a function of the target outage probability. 

Fig. 0] shows that the robust scheme achieves the target out¬ 
age probability for all channel amplitudes g. Consequently the 
same QoS can be obtained for all transmissions. In contrast, 
non-robust RA selects the transmission rate based Shannon 
capacity. CSI is assumed to be precisely available. Hence, the 
target outage probability is ignored by the algorithm. As a 
consequence, the requested QoS constraints are only achieved 
incidentally for small amplitudes. However, reducing the back¬ 
off factor a can increase the percentage of fulfilled QoS 
constraints. 

The throughput which results from the actually received 
transmission blocks is illustrated in Fig. |5] In the low-pout 
regime the guaranteed QoS comes with a reduction in through¬ 
put. For target outage probabilities between 0.1 and 1, the 
throughput decreases again. This effect results from a high 
packet error rate, which is forced by the robust algorithm 
in order to reach the target outage probability. However, this 
regime is of little practical relevance. 

IV. Robust Proportional Fair Scheduling with 
Imperfect CSI 

In this section, the results derived for robust rate adaptation 
are applied to the proportional fair (PF) scheduler. It is shown 
that the robust PF solution is achieved by just substituting 
the rate assigned with the one resulting from robust rate 
adaptation. Consequently, accounting for imperfect CSI does 
not change the scheduler itself but only its input. 

Furthermore, the optimal robust PF solution for delayed 
transmission acknowledgements is derived. In this case, the 
actual throughput during the last A transmissions is unknown 
to the scheduler as well. The optimal solution inherently 
estimates the unknown throughput by taking all possible cases 
into account. 

The previous section presented rate adaptation indepen¬ 
dently of the actual user and time slot, i.e., the respective 
indexes k and n are utilized again. In the following, we first 
present the model for the expected rate, we then discuss the 
case of immediate acknowledgements, and finally we present 


the robust PF scheduler for delayed transmission acknowledge¬ 
ments. 


A. Expected Rate 

Assuming UE k is served in time slot n, the expected rate 
is given by 

Rk[n] = E{i?fc[n]} = (1 - Pout,k)R[n], (16) 

which is the product of the probability of a successful 
transmission and the respective transmission rate. Since the 
experienced user rate is not precisely known to the transmit¬ 
ter, it refers to a random variable with mean value Rk[n]. 
Consequently, the actual experienced rate is given by 

Rk[n] = Rk[n] + 'll k[n], (17) 

which is the sum of the expected rate and a bias term 'I'fc[n]. 
The experienced rate is either the rate Rk[n] assigned for 
transmission or zero. Hence, the bias term results in 

\Rk[n]-Rk[n] if Rk[n] < Ck[n] 

\-4N o,he™„. 

By definition, the expectation of T'fe[n] is zero, since 


= {l-pout,k)iRk[n]-Rk[n])-poui,kRk[n\ 

= {l-pout,k)Rk[n]-Rk[n] = 0. 

With the given rate model the derivation of the PF scheduler 
can now be extended to the case of imperfect CSI, as it is 
stated in the following two sections. 


B. Immediate Acknowledgements 

In this section, we assume that the throughput of each 
user until the last time slot, n — 1, is precisely known to 
the transmitter. In practice, this would require immediate 
transmission acknowledgements. Although this case is rather 
artificial, it is an important step towards the more practical 
case of delayed acknowledgements. In the following, we use 
the indicator function Ik[n], which is equal to one if UE k is 
scheduled at time slot n and zero otherwise. 


Theorem 1. The robust PF scheduler selects the user with the 
largest value 


Uk[n] = K{uk[n]) 


Ik[N]Rk[N] 

Tk[N-l] 


( 20 ) 


Consequently, the PF scheduler itself does not need to be 
changed. Only the rate assigned for transmission is adapted 
according to the robust algorithm stated in Section \ni\ 


Proof The throughput of UE k experienced until time slot N 
is given as 


Tk[N] =jrJ:L,Ik[n]Rk[n] 

= ^ {Eu=i h[nmn] + 4[iV]i?,[A]) (21) 

= ^Tk[N - 1] + i4[iV]i?fe[iV] = 0. 
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Inserting (fTTT t into (ISTT t results in the expected throughput, 
given as 

E{Tfe[iV]} = ^^Tk[N - 1] + l4[iV]i?4Ar]. (22) 

Note, that (l22l) results from the zero mean random bias term, 
as it is given in (fT9l) . 


the largest value for 


z>fc[n] = E{z/fc[n]} 


Ik[N]Rk[N] 

Tfc[iV-l] 


(29) 


maximizes the expected utility increment. Note, that with (|29] | 
also the utility function M(E{T[iV]}) is maximized. □ 


The throughput values of all K UEs given in (l2Tl i are 
collected in vector 


T[N] = [T^[N],...,TK[N]f. (23) 

As dehned in m, the PF scheduler aims to maximize the 
utility function 


K 

u{T[N]) = '£^og{n[N]) 

k^l 

for an infinite number of time slots, expressed as 

K 

max lim log (T^ [TV]). 

k^l 

The gradient of (l24l i results in 

dTi[N] ’ ■ ■ ■ ’ OTkIN] _ 
{T,[N])-\...,{TK[N]r^' 


Vit(T[7V]) = 


(24) 


(25) 


(26) 


The second line comes from dlog{x)/dx = 1/x. Since the 
scheduler acts on time slot basis, the utility increment is of 
interest, which can be expressed by 

M(T[Af]) -M(T[Af- 1]) - 

- Vu iT[N - l])"^ {T[N] - T[7V - 1]) 

^ Tfc[7V] - rfe[7V - 1] 

Ip ^ Ik[N]Rk[N]-Tk[N-l] 

NTk[N-l] 

_^Ik[N][Rk[N] + ^k[N]) K 

NTk[N-l] N' 


The hrst and second line result from rearranging the difference 
quotient, while line three comes from inserting (|26] | and 12 ^ . 
With the dehnition of throughput in (12111 and the substitution 
of the actual transmission rate in (Ell the last two lines can 
be stated. 


Consequently, maximizing the utility increment dZTl i results 
in scheduling the UE with the largest value for 


C. Delayed Acknowledgements 

In this subsection, the robust PF scheduler derived in 
Section ITV-BI is extended to the more realistic case of delayed 
transmission acknowledgements. One of the main sources 
of CSI imperfections is a delayed feedback and backhaul 
transmission. Hence, also the confirmation about the success 
of previous transmission is not known instantaneously to the 
scheduler. 


Theorem 2. The robust PF scheduler for delayed transmission 
acknowledgements selects the user with the largest value 


Vk[n] 


Ik[N]Rk[N] 

fk[N] 


(30) 


where Tfc[-^] = 1/Sm=i denotes the throughput 

expected to be achieved at user k. In this regard the constant 
Wq = {N — A ~ 1)/A^ ■ Tfe[A — A — 1] is defined as well 
as the discrete random variable D = -A^k\j\Rk[j\, 

which comprises M = 2^ events, while v = 

Consequently, beside the rate assigned for the upcoming 
transmission also the value for the throughput which is as¬ 
sumed to be achieved by user k so far needs to be adapted. 


Proof. The throughput of UE k achieved up to time slot 
A is a random variable. First, it is known that UE k has 
been scheduled for v = '^n=N-A^k[n] transmissions within 
the period of interest. Hence, there are M = 2^ possible 
combinations of successful and non-successful transmissions. 
Each combination has a certain probability, resulting from the 
probability of outage for each transmission. The throughput 
equation (l2T1 i can be rewritten to 


Tm = mRk[i\+ 

+ Ef=”T-A 4b1i?fc[j] + 

= - A - 1] + 

+w Ef=T-A h[j]Rk[j\ + jrHN]Rk[N]. 

While Wq = {N — A — 1)/A • [A — A — 1] is known to the 
scheduler, the remaining part of the sum in (ISTT i is not, but can 
be split up according to (flTl i. Therefore, the utility increment 
as stated in (|27] | is given as 


z/fc[n] 


4[A](j?fc[A] + T'fc[A]) 
Tfc[A-l] 


(28) 


However, 'I'fc[A] is unknown to the scheduler and the ex¬ 
pected utility increment E{/(T[A]) — / (T[A — 1])} need 
to be maximized. As already stated in ( fT9] l the expectation 
EltkfclA]} is equal to zero. Hence, scheduling the UE with 


E{w(T[A])-m(T[A- 1])} = 

A-1 A 4[A]Afc[A] K (32) 
~ N ^ wqN + D a’ 

k—1 

where D = YI^Jn-a A is a discrete random variable 

with M = 2'^ events, each of which denoted as dm, Vm. 
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Consequently, maximizing the expected utility increment lead 
to scheduling the user with the largest value 


[n] = E 


wqN + D ) 


= E 

= E 


h[N]Rk[N] Jfc[jV]^fc[7V] ] 

WqN + D WqN + D j 

Ik[N]R^) 

WqN + D J 


(33) 


The third line in ([33]) is obtained, since the expectation w.r.t. 
the discrete random variable can be split up into a finite sum, 
where each summand is equal to zero, due to E{T'fc[n]} = 0. 
The expectation in the third line of (i3?t can also be calculated 
by summing over terms for the discrete events, weighted by 
the probability of occurrence, to 


M 

4[n] = Ik[N]Rk[N] Y, 

m—1 


1P{D = dm} 
WoN + dm 


(34) 


The expected throughput is given as Tfc [iV] = 
^ ^ immediate feedback, (i34l i 

also maximized the utility function u(E{T[7V]) for the 
delayed feedback case. □ 


V. Simulation Results 

In this section, we illustrate the efficiency of the derived 
scheduling schemes. We evaluate the performance based on 
Monte-Carlo simulations considering a simple scenario, where 
a single BS serves two users on orthogonal resources (one 
in each time slot). Note that two users refer to the simplest 
scenario although the basic effects are captured. The users 
are randomly placed within a circular area around the BS 
with radius 250 m which corresponds to an inter-site distance 
of 500 m. The respective mean channel gain results from 
Afc = /3dfc“, where a = 3.5 and /3 = is chosen 

according to 3GPP urban macro scenario ll^ . The simulations 
are done for 10,000 user drops, where for each drop 100 
subsequent channel realizations are generated. The receive 
SNR at the border of the serving area (at 250 m distance 
from the BS) is assumed to be either 5 dB or 10 dB. The 
impairment of CSI is generated by assuming a normalized 
coherence time of Tq = 10 transmission blocks. In an LTE 
system with a transmission block duration of 1 ms and a carrier 
frequency of 2.6 GHz, this refers to user mobilities of 10 km/h. 
Furthermore, we assume Np = 8 pilot signals are included 
per block and no quantization error is assumed (Q = oo). 
The channel uncertainty (MSE between actual channel and 
its estimate) as a function of the feedback/backhaul delay is 
illustrated in Fig. |6] Fig. [7] shows the portion of fulfilled target 
outage constraints (max. variance of 10 %). If the proposed 
robust rate allocation (RA) is applied, the outage constraints 
are always satisfied. However, the non-robust RA with and 
without backoff shows a significant gap. The slight increase 
for large delays is caused by channel estimation/prediction, 
which inherently assesses a higher channel uncertainty with 
smaller CSI amplitudes and reduces the outage probability. 



Fig. 6. Channel uncertainty (MSE between actual channel and its estimate) as 
a function of the normalized feedback/backhaul delay for two different SNR 
values. 
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Fig. 7. Portion of fulfilled target outage constraints as a function of the 
normalized delay. 


For each channel realization, at first rate adaptation is per¬ 
formed for each user and then one of the users is selected for 
transmission, according to the respective scheduling scheme. 
In the following, we compare proportional fair (PF) scheduling 
with perfect CSI (P-CSI) and imperfect CSI (TCSI), while for 
the latter case non-robust scheduling as well as the proposed 
robust scheduling is applied. The non-robust scheduler expects 
the available CSI to be perfect. Again, a back-off factor of 
a = 0.95 and a = 1 is used for simulations. 

For robust scheduling the two variants proposed in Section 
lIV-Bl and Section lTV-Cl are distinguished. The algorithm, which 
assumes immediate acknowledgements refers to Alg. 1. In this 
case, the scheduler just uses the acknowledgements available. 
For Alg. 2 the algorithm for delay feedback is applied, where 
an inherent estimation of the missing transmission acknowl¬ 
edgements is performed. 

The PF metric as it is defined in (l24l i as a function of the 





















































Fig. 8. Proportional fair (PF) metric as a function of the normalized delay 
for K = 2 users and SNR = 5 dB. 



Fig. 10. Resulting outage probability as a function of the normalized delay 
for Ff = 2 users and SNR = 5 dB. 
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Fig. 9. Propoifional fair (PF) metric as a function of the normalized delay 
for ii' = 2 users and SNR = 10 dB. 


Fig. 11. Mean user throughput as a function of the normalized delay for 
K = 2 users and SNR = 5 dB. 


normalized delay A is illustrated in Fig. Note that the 
PF utility metric u(T[7V]) = defined 

as the sum of logarithmic throughput values. The highest 
values are achieved for perfect CSI (black solid line), which 
is independent of the delay. Performing the same algorithm 
with imperfect CSI (without backing off) results in the green 
solid line. This scheme basically suffers from high outages 
(see Fig. [loll. Increasing the delay can even lead to marginal 
improvements in the PF metric. This is again a result from 
channel estimation/prediction as discussed before. Note that 
even for a delay of A = 0 blocks, the CSI is imperfect 
due to noisy pilot reception. Employing a back-off factor of 
a = 0.95 (green dashed line with ’o’ markers) causes a much 
higher percentage of successful transmissions, especially for 
small delays. Therefore, the PF metric significantly increases. 
However, even if the outage probability is very low in this 
region, it is not constant for all CSI amplitudes (see Fig. 
Si. Hence, the QoS constraints cannot be achieved at each 


transmission. 

The proposed robust PF algorithms (blue lines) outperform 
the non robust solutions for the SNR = 5 dB. However, 
the performance advantage disappears by increasing the SNR 
(compare Fig. [8] and Fig. S. For an SNR of 10 dB, small 
performance degradations can be observed for low (< 4 
blocks) and high latency (> 10 blocks). However, as can be 
seen from Fig. (TO] the target outage probability of pout = 0.1 
is achieved over the whole latency range and is particularly 
ensured at each transmission (see Fig. H. Employing Alg. 2 
(blue dashed line with ’o’ markers) results in small gains of 
the PE metric compared to Alg. 1 (blue solid line). This gain 
increases with A. 

In addition, the mean user throughput is illustrated as a 
function of the normalized delay (see Eig. fTTI i. In general, 
the throughput decreases with the delay. Small differences 
between the illustrated scheme can be observed, especially 
at delays between 5 and 15 blocks. Eor large delays, the non- 
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robust scheme wth back-off factor performs equally to the 
robust ones, while only the robust algorithm ensures the QoS 
constrains. In addition, the performance gain by employing 
Alg. 2 instead of Alg. 1 is rather marginal but increases with 
the feedback delay. 


VI. Conclusions 

In this work, a robust rate adaptation scheme together with 
a robust proportional fair (PF) scheduling algorithm have been 
presented, where both take into account that the available CSI 
is impaired. The algorithms aim to achieve a fixed outage 
probability at each transmission (independent of the channel 
amplitude), as it is of interest for delay critical applications. 
The robust PF algorithm has been derived for immediate as 
well as for delayed transmission acknowledgements. For the 
first case, the robust PF scheduler only requires a change 
in the rate adaptation, while the scheduler itself can remain 
unchanged. If the success of former transmissions is not 
immediately known at the BS side, the scheduler inherently 
performs an estimation of the already achieved throughput per 
user. Performance advantages of the proposed schemes have 
been shown for a two user example scenario. It was shown 
that PF utility metric as well as throughput gains increase with 
latency. System level aspects for scenarios with more than two 
users will be addressed in our future work. 


Appendix A 

Algorithmic Solution for Robust Rate Adaptation 


The rate adaptation function depends on three system 
parameters, the SNR 7 , the channel uncertainty e as well 
as the target outage probability pout- The numerical search 
for R based on (fTol i with the given system parameters and 
estimated channel amplitude g is quite complex. Hence, a real¬ 
time calculation might be impractical. An alternative solution 
is the pre-calculation of the rate adaptation function for a 
certain channel amplitude resolution. In this case, the real-time 
adaptation refers to a look-up table search. Such approximative 
scheme refers to discrete rate adaptation as discussed, e.g., in 

Ga. 

In the following, the calculation of (fT4l i is discussed by 
taking numerical limitations into account, such as, floating 
point resolution. First, it can be observed that the argument of 
the modified Bessel function in (fT4li decreases with the error 
e. Hence, small error values can result in large arguments and 
lead to values above the numerical resolution. In contrast, the 
exponential function in (fT4l i increases with the error e and 
hence balances the product. In other words, the product of 
the two factors might fall into the numerical resolution, while 
each factor may not. 

In order to overcome this problem, the modified Bessel 
function can be approximated for large arguments. The loga¬ 
rithm of the Bessel function behaves close to linear for large 
arguments, i. e., ln( Jm(x)) = amX + hm can be utilized: 


Jjnix) 


Jmix) X ^ ^rn 

exp(a + b 

m ) x>X m 



Fig. 12. The resulting outage probability as a function of the targe outage 
probability. 


The constants can be obtained by selecting two points Xi and 
X 2 and calculating the steepness according to 




\n{Jm{X2))-ln{J^{Xi)) 

X2 — Xi 


(36) 


while the shift is given by 


b^=\n{J^{Xi))-a^Xi. (37) 

With the approximation in dTSl l. the outage probability given 
in (fT4l) results in 

00 

Pout — 1 ^ ^ djyil (38) 

m—0 

with the approximated summand 

dm = exp (co) Jm{c 3 b), (39) 

including the coefficients cq = ci + 026 ^ -f m(ln(c 4 ) — ln( 6 )). 
Cl = -g^/(Ae), C 2 = -l/(Ae), C 3 = 2g/(Ae), C 4 = g and 
b = -l/p. 

In case the modified Bessel function reaches the numerical 
limit, (|39] | can be written as 

dm — exp (co + CLyu csb -I- b 

m ) if csb > X 

m ■ (40) 


This expression allows to perform the calculation within a 
numerical stable range. However, the approximation includes 
inaccuracies which may lead to approximation errors. Inaccu¬ 
racies can be observed for small target outage probabilities. 
Consequently, the robust algorithm would assign higher rates 
and the resulting outage probability exceeds the target, as 
illustrated in Fig. [121 However, this approximation error can 
be avoided by forcing the rate-assignment function to be 
monotonically increasing. 

In order to evaluate dHl numerically, the sum over an 
infinite number of arguments need to be restricted to M^ax, 
resulting in 

Mmax 

Pout ~ Pout — 1 ^ ^ dm ■ 

m—O 


(35) 


(41) 
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Note that a small value of Mmax can cause further inaccuracies. 
However, simulative experiments showed that M^ax = 150 
lead to a sufficient accuracy. 
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